Pairing fluctuation effects on the single-particle spectra for the superconducting state 
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Single-particle spectra are calculated in the superconducting state for a fermionic system with an 
attractive interaction, as functions of temperature and coupling strength from weak to strong. The 
fermionic system is described by a single-particle self-energy that includes pairing-fluctuation effects 
in the superconducting state. The theory reduces to the ordinary BCS approximation in weak 
coupling and to the Bogoliubov approximation for the composite bosons in strong coupling. Several 
features of the single-particle spectral function are shown to compare favorably with experimental 
data for cuprate superconductors. 
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Information on the single-particle spectral function, 
that is obtained from ARPES 1 and tunneling data 2 for 
cuprate superconductors, can shed light on the charac- 
teristic features of the superconducting state as well as 
on its connection with the unconventional normal state 
above the critical temperature T c . 

Most prominent among these features are the continu- 
ous evolution of a broad pseudogap structure from above 
to below Tc, 1 ' 2 the emergence of a coherent peak below 
T c that combines with the pseudogap structure to yield 
a characteristic peak-dip-hump profile 1 (for which two 
distinct energy scales can be identified), and the pecu- 
liar dependence of the frequency position and weight of 
the coherent peak on temperature and doping. 1 Gener- 
ally speaking, features of the standard BCS theory are 
recovered for overdoped samples, while non-BCS behav- 
iors occur for optimally-doped and underdoped samples. 

The origin of the peak-dip-hump profile has especially 
been the subject of controversy, being attributed either to 
"extrinsic" effects like the bilayer splitting 3 or to "intrin- 
sic" effects. The latter can be identified over and above 
the extrinsic effects, and are believed to originate from 
strong (many-body) interactions in the system 4 . 

Quite independently from the microscopic origin of 
the fermionic attraction giving rise to superconductiv- 
ity, its strength is believed to be stronger in optimally- 
doped and underdoped than in overdoped samples , con- 
sistently with the above findings. This implies that pair- 
ing fluctuations should definitely be taken into account, 
irrespective of other effects (such as the bilayer splitting 
and/or additional many-body effects associated with spe- 
cific pairing mechanisms 6 ). 

In this paper, we assess the role of pairing fluctua- 
tions for the single-particle spectral function in the su- 
perconducting state on rather general grounds, by iden- 
tifying a single-particle self-energy that describes fluctu- 
ating Cooper pairs in weak coupling and non-condensed 
composite bosons in strong coupling. (The latter form as 
bound-fermion pairs in the strong-coupling limit of the 



fermionic attraction.) To this end, we consider fermions 
mutually interacting via an attractive contact potential 
in a 3D continuum, without taking into account lattice 
effects nor the explicit physical mechanism responsible 
for the attraction. Qualitative comparison with experi- 
mental data will thus rest on identifying the increasing 
coupling strength in this model with the increasing po- 
tential strength for decreasing doping level in the phase 
diagram of cuprate superconductors. Although the ef- 
fective model here considered is oversimplified for a full 
description of cuprates, the physical questions we are ad- 
dressing are sufficiently general that this minimal model 
will prove sufficient to capture the main experimental 
features. 

Our main results for the single-particle spectral func- 
tion (to be compared with the experimental findings) are 
summarized as follows: 

(i) A broad pseudogap structure and a coherent peak are 
simultaneously present in a wide range of coupling and 
temperature below T c , giving rise to a peak-dip- hump 
profile. Two distinct energy scales (the positions A pg of 
the broad pseudogap structure and A m of the coherent 
peak) can then be extracted from the spectra at given 
temperature and coupling, as seen experimentally 1 . 

(ii) At a given (notably, intermediate) coupling, BCS-like 
features coexist with non-BCS behaviors. The position 
A m and weight z of the coherent peak versus wave vec- 
tor follow a BCS-like dependence, as found experimen- 
tally 7,8 . At the same time, the weight z of the coherent 
peak has a strong temperature dependence, in agreement 
with experiments 9,10 but in contrast with BCS theory. 

(iii) At low temperature, the weight z of the coherent 
peak has a strong dependence on coupling, decreasing 
monotonically from weak to strong coupling (as seen ex- 
perimentally in cuprates for decreasing doping 9,10 ). At 
the same time, the position A m of the coherent peak in- 
creases monotonically with coupling 9 ' 10 . 

(iv) The positions A pg of the broad pseudogap structure 
at T c and A m of the coherent peak near zero temperature 
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cross each other as function of coupling about interme- 
diate coupling. This feature is also seen experimentally 
by intrinsic tunneling experiments at different dopings 11 
(although these data are subject to controversy 12 ). 

Pairing fluctuations are taken into account in our the- 
ory by considering, besides the off-diagonal BCS-like self- 
energy, the diagonal i-matrix self-energy suitably ex- 
tended to the superconducting state (we use the Nambu 
formalism throughout). The t- matrix self-energy has 
been widely used to include pairing fluctuations in the 
normal state 13 , and specifically to account for pseudogap 
features in the single-particle spectral function 14 . In our 
theory, the diagonal ^-matrix self-energy (that survives 
above T c ) will essentially be responsible for the pres- 
ence of the broad pseudogap structure below T c . The 
off-diagonal BCS-like self-energy will instead give rise to 
the simultaneous emergence of the coherent peak. In ad- 
dition, our theory recovers the Bogoliubov approximation 
for the composite bosons in strong coupling. 15 

The diagonal t-matrix self-energy reads 
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while the off-diagonal self-energy has the BCS-like form 
£12 (k) = —A in terms of the superconducting order pa- 
rameter A. The pairing-fluctuation propagator Tn in 
the broken-symmetry phase entering Eq. (1) is given by: 
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In the above expressions, q = (q, f2„), k = (k, u>i), and 
p = (p, u n ) (q, k, and p being wave vectors, f2„ a bosonic 
Matsubara frequency, and loi and u„ fermionic Matsub- 
ara frequencies), m is the fermionic mass, (3 the inverse 
temperature, and Gij = 1,2) have the form of the 
BCS single-particle Green's functions in Nambu notation. 

For given temperature and strength of the point- 
contact interaction, the chemical potential \x and the or- 
der parameter A are obtained by solving the coupled 
particle number and gap equations: 
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where rj = 0+ and E(k) = [£(k) 2 + A 2 ] 1 / 2 with £(k) = 
k 2 / (2m) — p. The scattering length ap of the (fermionic) 
two-body problem has been introduced to eliminate the 
ultraviolet divergences originating from the use of a 
point-contact interaction. The dimcnsionlcss parame- 
ter (kpap)^ 1 (where kp is the Fermi wave vector re- 
lated to the density n via kp = (3tt 2 ti) 1/ ' 3 ) character- 
izes the fermionic coupling strength and ranges formally 
from — oo to +oo. [Correspondingly, the Fermi energy 
Ef equals k F /(2m).] In practice, the crossover from the 
weak- to strong-coupling regimes occurs in the limited 
range —1 < (kpap)^ 1 < +1, which will be mostly ex- 
plored in what follows. The dressed normal (diagonal) 
Green's function Gn in Eq. (5) is obtained by matrix 
inversion of the Dyson's equation in Nambu formalism 
with the above self-energies and is given by: 
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Note that, while the number equation (5) contains 
the dressed normal (diagonal) Green's function, the 
gap equation (6) is obtained from the anomalous (off- 
diagonal) BCS Green's function and thus its form is not 
modified with respect to the BCS theory. This ensures 
that the pairing-fluctuation propagator (2) remains gap- 
less for all temperatures (below T c ) and couplings. [The 
value of T c is obtained from Eqs. (5) and (6) by setting 
A = identically] The numerical values of the chemi- 
cal potential \i and order parameter A at given tempera- 
ture and coupling differ, however, from those obtained by 
BCS theory due to the different structure of the number 
equation. It can be verified 15 from Eqs. (2)-(6) that the 
Bogoliubov approximation for the composite bosons is re- 
covered in the strong-coupling ([3/j, — > — oo and A <C |/x|) 
limit. This represents a notable achievement of our the- 
ory. 

The use of bare BCS single-particle Green's functions 
in the self-energy (1) further enables us to perform the 
analytic continuation to real frequency uo in a closed form, 
as to avoid numerical extrapolation procedures. To this 
end, and similarly to what was done in Ref. 14 for the 
normal phase, we express the pairing-fluctuation propa- 
gator Tn by its spectral representation: 
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Here ImFn(q, lu) is defined as the imaginary part of 
rn(q, iCl v — > u + it]). This quantity is, in turn, ob- 
tained from the expressions (2)- (4), with the replacement 
iQ, v — > cj + ii] made after the sum over the internal (Mat- 
subara) frequency has been performed. In this way, the 
imaginary part of the retarded self-energy is: 
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Im£f 1 (k,w) = -y ^{ u 2_ k i m r 11 (q,w + £;(q-k)) 

x [/(£(q-k)) + 6(w + ^(q-k))] 

+ ^_ k Imrn(q,w-£7(q-k)) 

x [/(-£7(q - k)) + 6( W - £7(q - k))]} . (9) 

Here, f(x) = [exp(#c) + l]" 1 and 6(x) = [exp(/3a;) - l]" 1 
are the Fermi and Bose distributions, while and u k 
are the usual BCS coherence factors. The real part of 
S{\ is then obtained via a Kramers-Kronig transform. 
The retarded single-particle Green's function Gf\ results 
from these functions, in terms of which the single-particle 
spectral function A(k.,uj) = — -ImG^(k, uj) of interest is 
eventually obtained as a function of lo for any given k. 16 
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FIG. 1. Single-particle spectral function vs frequency (in 
units of sf) at different temperatures, for |k| = Ay and 
(fcFaF) -1 = -0.45. 

Figure 1 shows the evolution with temperature of the 
single-particle spectral function for an intermediate cou- 
pling ((kpap)^ 1 = —0.45). [For this coupling, the ra- 
tio of the pair-breaking temperature T* to the critical 
temperature T c is about 1.25, as obtained in Ref. 14.] 
Focusing specifically on the features at negative ui, note 
how the (sharp) coherent peak grows from the (broad) 
pseudogap structure already present at T c . The coher- 
ent peak becomes sharper upon lowering the tempera- 
ture and gains weigth at the expenses of the pseudogap 
structure, giving rise to a characteristic peak-dip-hump 
profile. The two features coexist over a wide range of 
temperature. 

Figure 2 shows the wave-vector dependence of the 
single-particle spectral function at the temperature T = 
0.7T C for the same coupling of Fig. 1. Following the evolu- 
tion of the coherent peak across the underlying Fermi sur- 
face, one identifies the characteristic particle-hole mix- 
ing of the BCS theory, with a reflection of both particle 
and hole bands accompanied by a transfer of the spectral 
weight from negative to positive frequencies. In our the- 
ory, such BCS-like features coexist with non-BCS behav- 
iors, namely, the occurrence of the pseudogap structure 



and the strong temperature dependence of the spectral 
weight of the coherent peak (see Fig. 4 below) . 
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FIG. 2. Single-particle spectral function vs frequency (in 
units of Sf) at different wave vectors about fc p / , for T = 0.7T C 
and (/cf^f) -1 = —0.45. 



0.5 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Mf) 

FIG. 3. (a) Weight z and (b) position A m (full line) of 
the coherent peak at negative frequencies vs coupling for 
T = 0.1T C . In (b) the pseudogap A pg for T — T c (dashed 
line) is also shown. 

Figures 3(a) and 3(b) report, respectively, the weight 
z and position A m of the coherent peak of the single- 
particle spectral function at negative frequencies vs cou- 
pling for a low temperature (T = 0.1T C ). The peak 
weight z saturates at the BCS value 0.5 in weak cou- 
pling, decreases markedly across the crossover region 
— 1 < (kpap)^ 1 < +1, and becomes negligible in strong 
coupling. At the same time, the peak position A m in- 
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creases monotonically across the crossover region. In 
Fig. 3(b) we also report the pseudogap A pg (dashed 
line), as identified from the position of the maximum of 
the spectral function at T c . While the qualitative trend 
of A m and A pg vs coupling is similar, the two curves 
cross each other at about the intermediate-coupling value 
{UfCLf) 1 = —0.45. In addition, we have verified that 
A m about coincides with the value of the order parame- 
ter A in the weak-to-intermediate coupling region. 
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FIG. 4. Temperature dependence of the weight z of 
the coherent peak at negative frequencies (full line) for 
(fefCiF) -1 = —0.45. The superfluid density p s (dashed line) 
is also shown for comparison. 

Figure 4 reports the dependence of the weight z on 
temperature for the coupling value (fc^ftF) 1 = —0.45 
(full line). Note the strong temperature dependence of 
this quantity, which vanishes at T c . This contrasts the 
BCS behavior, whereby the weight of the coherent peak 
for negative frequencies would equal 0.5 irrespective of 
temperature. The temperature dependence of the su- 
perfluid density p s (calculated according to Rcf. 17) is 
also shown in the figure (dashed line). The resemblance 
between these two quantities has been noted experimen- 
tally for nearly-optimally-doped cuprates 9 ' 10 . We have, 
however, verified that this resemblance docs not occur 
for other values of the coupling, both on the weak- and 
strong-coupling sides of the crossover. On the weak- 
coupling side of the crossover, in particular, our finding 
is confirmed by the BCS theory, whereby z = 0.5 irre- 
spective of temperature while p s decreases monotonically 
from T = to T c . For this reason, no universal corre- 
spondence between the temperature dependence of z and 
p s should be expected on physical grounds. 

The results shown in the above figures refer mostly to 
the weak-to-intermediate coupling side of the crossover 
region. Our conclusions have, however, been drawn from 
a sistematic study of the whole crossover region from 
weak to strong coupling, of which Fig. 3 is an example. 

All the qualitative features extracted from the above 
figures compare favorably with the experimental data on 



cuprates, as anticipated in the points (i)-(iv). This qual- 
itative comparison rests on the assumed correspondence 
between the increasing of the coupling strength in the 
present theory and the decreasing doping level in the 
phase diagram of cupratc superconductors. 

In conclusion, we have shown that pairing fluctuations 
can (at least qualitatively) account for several nontrivial 
features of single-particle spectra in the superconduct- 
ing state. Our results specifically demonstrate that the 
experimental finding of two distinct features (pseudogap 
structure and coherent peak) in the single-particle spec- 
tral function is fully consistent with the occurrence of 
strong pairing fluctuations in cuprate superconductors. 
Competition of two distinct order parameters is there- 
fore not required to account for the occurrence of two 
different energy scales in the experimental data. 
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